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Abstract 

Considering the hydrodynamical limit of some interacting particle sys- 
tems leads to hyperbolic differential equation for the conserved quantities, 
e.g. the inviscid Burgers equation for the simple exclusion process. The 
physical solutions of these partial differential equations develop discon- 
tinuities, called shocks. The microscopic structure of these shocks is of 
much interest and far from being well understood. We introduce a domain 
growth model in which we find a stationary (in time) product measure for 
the model, as seen from a defect tracer or second class particle, traveling 
with the shock. We also show that under some natural assumptions valid 
for a wider class of domain growth models, no other model has stationary 
product measure as seen from the moving defect tracer. 

Key-words: second class particle; shock solution. 

Introduction 

The hydrodynamical limit of the nearest neighbor asymmetric simple exclusion 
model leads to the inviscid Burgers equation 

du 1 du 2 
~dt + 2~dx~ = 

which is a special case of the one-component hyperbolic conservation law 

where u i— > J(u) is a smooth, typically convex function. (By changing x to —x, 
concave J-s can be transformed to convex ones.) This equation has a shock 
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(weak) solution starting with initial data 



u(0, x) = 
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with uieft > u r ight ■ The stable weak solution is of the form 



u(t, x) = 



( 
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x > st 



where the speed s of the traveling shock is determined by the Rankinc-Hugoniot 



see e.g. [13]. This is what we see on a macroscopic scale. The microscopic 
structure (i.e. on the level of particles) of the shock is of great interest. It has 
been considered in the context of the asymmetric simple exclusion process, and 
rather complicated microscopic structures have been found [2] [3] [4] [5] [6] . In 
the more general context of attractive particle systems the microscopic structure 
of the shock was investigated by [11]. 

In the present note we consider a class of one-dimensional domain growth 
models, parametrised by a jump rate function, r : Z — > M. In a special case of 
the rate function we show that the shock, as seen from a defect tracer (second 
class particle) has stationary (in time) distribution of product structure which 
we identify. We also show that this is a peculiarity of the case considered, no 
other model in the wide class of these models has this property. The structure 
of the paper is the following: 

In section 2 we define the class of models considered and determine the station- 
ary distributions for them. 

We describe the hydrodynamic limit of these models and calculate the speed of 
the shocks using Rankine-Hugoniot formula (2) in section 3. 
In section 4 we introduce the defect tracer in our models. Via Rankinc-Hugoniot 
formula, we also give an indication on the fact that, in general, shock solutions 
are closely related to measures stationary as seen from the defect tracer. 
The last section contains our main result on the product structure of such a 
stationary distribution as seen from the defect tracer. This gives an explicit 
description of the microscopic shape of some types of shock solutions. 
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■/("right) - J(Mlcft) 
Wright _ Ui e ft 
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2. The bricklayers' model 



2.1 Infinitesimal generator 



We consider the phase space 



= {uj_= (uji) ie z : Ui £ Z} = Z' 



■z 



For each pair of neighboring sites i and i + 1 of Z, we can imagine a column 
built of bricks, above the edge (i, i + 1). The height of this column is denoted 
by hi. If £ O for a fixed time tel then Ui(t) = /ij-i(i) - G Z is the 
negative discrete gradient of the height of the "wall" . The growth of a column 
is described by Poisson processes. A brick can be added to a column: 

{uji, — ► (ui - 1, + 1) with rate r(u>i) + r(-oj i+1 ) . 

See fig. 1 for some possible instantaneous changes. The process can be repre- 
sented by bricklayers standing at each site i, laying a brick on the column on 
their right with rate r(u>i) and laying a brick to their left with rate r(— u>i). This 
interpretation gives reason to call these model bricklayers' model. For small e 
the conditional expectation of the growth of the column between i and i + 1 in 
the time interval [t, t + e] is {r(wi(t)) + r(—Ld i+ i(t))} ■ e + o(e). Note that the 
process has a left-right mirror symmetry, i.e. the rate of a column's growth is 
the same as if looking at the reflected configuration. We want the dynamics to 
smoothen our interface, that is why we assume monotonicity of the rate function 
r, which means that a column grows more rapidly if it has a higher neighbor 
on the right or on the left. In later sections we shall impose another restrictive 
condition on r, see (5). 

At time t, the interface mentioned before is described by to(t). Let <p : 
£1 — > K be a bounded cylinder function i.e. <p depends on a finite number of 
values of Wj. The growth of this interface is a Markov process, with the formal 
infinitesimal generator L: 



{L<p)(tj) = [r(wi) + r(-w i+ i)] • [<p(. . . , m - 1, uj 1+1 + 1, . . . ) - (p(u)} 



Note that for each index i, Wj can also be negative hence direct particle inter- 
pretation fails, see the remark after formula (4). 

When constructing the process rigorously, problems may arise due to the 
unbounded growth rates. The system being one-component and attractive, we 
assume that existence of dynamics on a set of tempered configurations f2 (i.e. 
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r(u) t ) +r(-ui t+1 ) 



i ffl 



Figure 1: A possible move 
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configurations obeying some restrictive growth conditions) can be established 
by applying methods initiated by Liggett and Andjcl [1] [7]. Technically we 
assume that is of full measure w.r.t. the canonical Gibbs measures defined in 
2.2. We do not deal with this question in the present paper. 

The exponential bricklayers' model 

A special case of the models is the exponential bricklayers' model (EBL), where 
for z E Z 

(3) r(z)=e"5^ 



with a positive real parameter /3. 



2.2 Translation invariant stationary product measures 

In this subsection we show a natural way to construct a stationary translation 
invariant product measure for our models. By chapter one of [9], a measure /i 
is stationary iff for any bounded cylinder function <p, 

E(Lip)(u) = 

is satisfied for a process distributed according to \i. We assume /i to be a product 
measure with marginals 

H{z) = fi{uJ. : u>i = z} 
for zeZ. By changing variables and using product structure of //, 



E(L<p)(u) = 

= E 5Z{[ r ^) + r (- w i+l)] • [<P(- •■ , <*>i - 1, + 1, . . .) - <p(w)] j = 



iez, 



r(wi+l) 



+ 1) /J,(u) i+1 - 1) fj,(wi + 1) - 1) 



fj,(ui) M(^i+i) 



-+r(-w i+ i + l)- 



/x(Wi) M(<*>»+i) 



- r(wi) - r(-w i+ i) 



This expression becomes zero if we make the sum telescopic on the cylinder set 
supporting if. Hence stationarity of jj, is assured by assuming 



r(z) 



H(z - 1) 



fi(z) 



— and r(—z) ■ — , 

1) V ; (i(z + l) 
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to be constants. As a consequence, we obtain the condition 



(4) r(z) ■ r(—z + 1) = constant . 
There are two essentially different choices. 

r(z) ■ r(-z + 1) = 

defines models of zero range types, we do not consider this possibility here. The 
other choice is choosing the right-hand side of (4) to be a positive constant. In 
this case, by rescaling time, we can turn this constant to be one without loss of 
generality: 

(5) r(z) -r(-z + l) = 1 . 

Rates (3) of the EBL model satisfy this condition. 
For n E N, we define 

n 

r{n)\ := J| r(y) 

!/=l 

with the convention that the empty product has value one. Let 

9 : = log ( liminf (r(n)!) 1 ^™ ) = lim log(r(n)) , 

which is strictly positive by (5) and by monotonicity of r, and can even be 
infinite. With a generic real parameter 9 e (— 9, 9), we define 

OO n — 

Z{6):= £ 



Z— — OQ 



r{\z\)\ 

and the product measure ji^ with marginals 



(6) ^ '(*):= ' 



Z{9) r(\z\)\ 
which has the property 

r( - z) '^TT) =c 
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thus it is stationary. We call these measures canonical Gibbs-measures. 

For the special case of the EBL model, for 9 E (—00, 00), we obtain the 
discrete normal distribution 

(8) M (8) (*) = — 55 = ~ - 

e -wZ(6) Z(f3, m) 

with the notation m := 9/(3. 

3. Hydrodynamical limit 

Being attractive due to monotonicity of r, we can take the hydrodynamical limit 
of a bricklayers' model by 

(9) u(t, x) :=Ew x/e (i/e) • 

Then via formal computations we obtain differential equation (1) 

du dJ(u) 
1 — = 

dt dx 

as e — > 0, where J(u) is defined as follows. The function u{9) = E^(uj) of 9 is 
strictly increasing since the derivative 



^ = (e<«V)-(e<*>m) 9 ) 



is positive (—9 < 9 < 9). Let 9{u) be the inverse function. The quantity 
(r(cj) + r(—u)) depends on 9, and J is defined by 

(10) J{u) := E (9(M)) (r(w) + r(-w)) = 2 cosh(#(u)) . 

Proposition 3.1. There exist 9\ < < 62 numbers such that J(u) defined 
above is convex on the interval (u(9i), 14(^2))- 

Proof. With the notations 

«W : =-dT and u W : ="d^- 

and by computing derivatives of inverse functions, we obtain from (10) 

d 2 J cosh(6>) 1_ sinh(fl) m"(0) 

d^2°( u W)- 2 (u'W) 2 2~ (u'(0)) 3 ■ 
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The positivity of the left-hand side is assured in case 8 = 0, and is equivalent 
to the condition 



(11) 



< ctanh(0) , if > 

u'{0) 
u"{6) 



u'(0) 



> ctanh(6>) , if < 



by positivity of u'(0). The function i— ► u(0) is analytic in (—6, 0), u'(0) is 
strictly positive, hence the left-hand side of (11) is bounded on the interval 
(-0*, 8*) for each < 0* < 0. Due to the unbounded behavior of ctanh(#) on 
any interval containing zero, there exist 8\ < < 6* 2 for which (11) and hence 
convexity of J(u) is satisfied. □ 

Using definitions (10) and (9) in Rankine-Hugoniot formula (2), the speed 
of the shock can now be written as 

_ E( 8 ( U " E ")) (r(w) + r(-Lo)) - EW"><*» (r(Lu) + r(-uj)) 
^ ' S ~ E( 9 ("" g ht))( w ) - E( 9 («^t))(w) ' 

4. The defect tracer 
4.1 Coupling the models 

Let lj + (0) and w~(0) be two elements of £1 At time t = wc start with a 
configuration where these two realizations differ at only one site: 

w^(0)=wr(0) if i^O and lu+ (0) = lu {0) + I . 

One possible representation can be imagined by two walls. At time 0, the walls 
are the same on the right side of position 0, and every column of the wall + is 
higher by one brick than column of wall ~ on the left side of zero. We want the 
two processes to grow together in such a way, that the difference between them 
remains "one step" at any time t > 0: 

(Vt > 0) (3i Q(t) G Z) : w+(t) = w r(t) if i ^ Q(t) and 

W Q(t)(*) =W Q(t)(*) + 1 • 

We shall call this difference between the two models defect tracer, and Q{t) is 
its position at time t. We show the coupling which preserves the only one defect 
tracer while both lu~ and cj + evolves as usual. This coupling for the simple 
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with rate 




*0-it 






Q has . . . 




• 








decreased 


r ( w Q-l) + r (" w Q) 


• 


• 
















• 


increased 


K w q)+K- w q+i) 






• 


• 





Table 1: The coupling rules 
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exclusion model is described (with particle notations) in [8] and [9]. Let our 
defect tracer be posed at point Q (i.e. ujq = ojq + 1; ojf = w t ~ if i ^ Q), 
and let hf | (or |) mean that the column of lj + (or the column of lo~, 
respectively) between the points i and i + l has grown by one brick. Then 
the growing rule for the columns hg_ 1 and Kq is shown in table 1. Every line 
of that table represents a possible step with rate written on the first column. 
These rates are non-negative due to monotonicity of r. For each column of this 
table, summing the rates corresponding to the possible steps assures us that 
columns of each lj + and u/~ evolve as usual in the neighborhood of Q. For 
i ^ Q — 1 or Q, h~l and ftr increases at the same time with the original rate 
+ r(-uj7 +1 ) = r {uj+) + r(-w+ 1 ). 
How does an observer following the defect tracer see the surface? We intro- 
duce the drifted form Tk w of an w € ft as follows. Let k G Z, then Ttw £ 
and 

(T fe w)i : = Wi_ fe . 

From now on, we denote by w(t) the w~ (i) process as seen from the position 
Q(t) of the defect tracer, i.e. Wj : = According to the coupling rules, we 

can write the infinitesimal generator for w: 

(13) (L^'VXw) = 

i#-l,0 

+ K w -i) + r (-^o - 1)] • [<p(- ■ ■ , w_i - 1, w + 1, . . . ) - <p(u)] + 
+ [r{u ) + r(-wi)] • [<p(. . . , w - 1, wi + 1, . . . ) - ¥>(w)] + 
+ I r (-^o) - r(-w - 1)] • [<p(ri(. . . , w_i - 1, w + 1, •••))- ¥>(w)] + 

+ [r(w + 1) - r(w )] • [<p{t-! w) - <p(w)] . 

4.2 The speed of the defect tracer 

The main problem of this note is to find a stationary measure for the process 
as seen from the defect tracer, i.e. to find a measure ^ d t ), for which 

E(L( d - t 'V)(w) = 

is satisfied. Before giving a partial answer to this question, we give an early 
indication on the correspondence to shocks of such a measure fi^ d,t '\ 

Let a < — 1 (and b > 1) be sites far on the left side (and far on the right 



10 



side, respectively) of the defect tracer. We choose the function 

b 

k— a 

in (13) to obtain 

(L^'VXw) - [r(w„_i) + r(-w„)] - [r(w 6 ) + r(-w&+i)] + 
+ [r(-cuo) - r(-cj - 1)] • (w a _i - w&) + [r(w + 1) - r(ui Q )] ■ (uj b+1 - u> a ) ■ 

Let us assume that a measure * ) is stationary for !,(<*•*•). Let us also assume 
that as I — > ±oo, the random variable w; becomes asymptotically independent 
of wo, u)\, and the distribution of loi converges weakly. Then we have 

(14) = E(L^<p)(u) = E [r(w„) + r(-w„)] - E [r(w 6 ) + r(-w 6 )] + 
+ E [r(-w ) - r(-w - 1)] • E(w a - w 6 ) + E [r{u + 1) - r(w )] • E(w 6 - uj a ) + 

+ H(a,b) , 

where the error function H(a, b) tends to zero if a — ► — oo and — > oo. (For 
the product measure ^t( d t ) we find in the next section, H(a, b) = for any 
a < — 1, 6 > 1.) According to the rules of the coupling, and assuming also 
ergodicity of the process as seen from the view of the defect tracer, we have the 
law of large numbers 

v : = lim Qlp- = E {[r{uj + 1) - r(w )] - [r(-wo) - r(-w - 1)]} a.s. 

for the speed of the defect tracer. Hence we conclude from (14) that 

E [r(u b ) + r(-u b )] - E [r(uj a ) + r(-w a )] 

a^-oo,b^oo E(W(,) — E(w a ) 

in case E(w ) 7^ E(w b ) and their limits are not equal i.e. the slope of the surface 
is different far on the two sides. This formula is the same as (12), which we 
obtained for the speed of the shock using the Rankine-Hugoniot formula. This 
shows that a measure /x( d t ) with different asymptotics on the left and on the 
right can be identified as the microscopic structure of a shock solution of (1). 
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5. Stationary measures as seen from the defect 
tracer 



In this section we find a stationary product measure satisfying (7) for the defect 
tracer of the EBL model. We also show that this kind of measure only exists 
for the EBL model. 

Intuitively one expects that far from the defect tracer a stationary measure 
behaves like the canonical Gibbs-measure since the defect tracer is only a 
local "error" for the evolution of the process. The canonical measure has one 
parameter 6, but it is not necessary that in this case parameter 6*i e ft far on the 
left side would be equal to the parameter 6> r i g ht far on the right side. Let 

: = {6i : i e Z} 

be a sequence of parameters. Then it seems to be reasonable to assume that 
the product measure /j,(-> with marginals 



Z(8 t ) r(\z\)\ 



is stationary for L d t - (13) (i e Z). This measure only differs from the canonical 
^ e "> (6) in that the parameter of its one-dimensional marginals depends on the 
position. The question is whether there are any choices of 8 for /j,^ to be 
stationary. 

Theorem 5.1. For a bricklayers' model, if r is not the constant function, then 
the measure fi@> described above is stationary for L( dd ^ if and only if r is the 
rate of an EBL model with any parameter [3 > 0, and for the 8 parameters of 
(£) 



ft 



Ouft if i < -1 , 

Omght ■ = Oleft ~P if % > 



is satisfied with an arbitrary real number 8i e ft- 
Proof. Stationarity means 

E(®(L d - t - V )(w) = , 
after some changes of variables, by straightforward computations we obtain from 
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(13) 



(15) = E^j{i4 + B + C' + £>} ¥ >(w)j , 

where 

+ 1) fj, i+1 (u} i+ i - 1) 



, Hi(ui) fj, i+1 (u i+1 ) 



A = Yl [[r(wi + 1) + r(-o; i+ i + 1)] 

- [r(wi) + r(-w i+ i)] , 

B = r(w_i + 1) + r(-w ) ? : — -—^ 

- r(u;_i) - r(-w ) - r(">o + 1) + r-(wo) > 

C = [r(-m + 1) - r(-a*)] • t \ } 1} • II ^TY ■ 



£) = [r(w_i + l)-r(w_i)]- J] 



We eliminate the expressions \i k for all k G Z with the use of (5) and (7) 



to obtain 
i4 = 



53 [e 9 --^ 1 r(wi+i) + e 9 '-^ 1 r(-wO - r(w<) - r(-w i+ i) 

r-(wo) 



S = c e - 1 - yo r( Wo )+c y - 1 - eo r(- W _ 1 ). 

r(w + 1) 

- r(u)-i) - r(-ujo) - r(uj + 1) + r(u ) , 

c = rM(i--^)^n 9(,, - H,H #i - 

V r(wi + l)y 

D = [r(,- 1 +l)-r(,_ 1 )]Ile^- ( >A . 

jez ^ 

For a < — 1, 6 > 1 fixed, let us consider bounded cylinder functions (p, which 
depend on the variables ix> a , w +ij • • • > By stationarity of /i® , (15) is satis- 
fied for all of them. Hence it is necessary that A + B + C + D does not depend 
on the variables u a , co a +i, ■ ■ ■ > and its mean is zero according to Only 
C, D, and the second term in B are the terms which can contain product of 
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functions of different variables ujk- Each of them is positive by monotonicity of 
r. Thus it follows that each of these three terms must not depend on more than 
one variable. This implies 



r(z + l) 



= constant = r(0) 



due to the form of the second term in B. The value r(0) 2 of this constant is a 
consequence of (5). Thus we conclude that r is necessarily exponential, the rates 
are that of the EBL model (3). C and D can also contain at most one variable, 
hence we obtain 6k = Q-i for k < — 1 and 9k = 9o for k > 0. This means that 
we have at most two kinds of marginals of one on the left-hand side of the 
defect tracer and an other one on its right-hand side. In (15), computing the 
expectation value of p times the terms of A, summed up for i < a — 2 and for 
i > 6+1, gives zero. The reason for this is that the variables in these terms are 
independent of the variables which ip depends on. For the rest of the indices, 
note that we have a telescopic sum for A. Due to this and using the rates of 
the EBL model, we can simplify our expressions to 

A = r(-u a -i) - r(w a -i) + r(uj b+ i) - r(-u> b+ i) + 

+ r(w-i) - r(-uj-i) + r(-wo) - r(u>o) , 

B = e 9 - 1 - 0o r(wo) + e £ '- 1 - eo -' 3 r(-cc ) _ 1 ) - 

- r(w_i) - r(-w ) - e r(u o ) +r(ui ) , 

C = (l- e ^)r(-co)e^-^^^\ , 

D = (e^-lJrMe^-O--.^ . 

It is easy to check that simply choosing 0_i = 9 does not eliminate the variables 
u>-i, luq from A + B + C + D. Hence this can not be a solution for 9 to make 
equation (15) be satisfied for all ip. This means that the marginals on the left- 
hand side of the defect tracer are different from those on the right-hand side. 
When taking expectation value in (15), this leads to having constant times <p 
from the terms containing uj a -i and uib+i in the first part of the expression of 
A. In order to make (15) be satisfied for all ip, it is necessary that we obtain 
other constants to have zero together with. They can only come from C and D. 
Thus we conclude 

(16) e e -i- 6 ° = 



14 



with the use of the form (3) of r. In view of (8), we have the measures 



P 2 V p J p — 2 

(17) m(z) = - 



e""2(«i) 



Z(/3, mi) 



with mi := #i//3. We know that the normalization Z((3, m) in the right-hand 
side of (17) is periodic in the parameter m with period one, which tells us 



Z(g_i) _ Z(fl + g) _ e ^ g + 1) 

Z(8 ) - Z(6 ) 6 ) 



Using this result together with (16) and with the property E'-V(±Wi) = e ±9i , 
we see that (15) is satisfied, which completes the proof. □ 

The form of the measure described in this theorem shows that the discrete 
normal distribution of u>i, i < — 1 is shifted by +1 compared to the distribution 
of Wj, j > 0. This gives us the picture of a (random) valley with the (randomly) 
moving defect tracer in its center. Since the position of the defect tracer is not 
deterministic, we do not see the sharp change between the distribution of the 
two sides of this valley, if looking the model from outside. 
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